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ON THE BI-HAMILTONIAN THEORY FOR THE HARRY DYM 

EQUATION 

MARCO PEDRONI, VINCENZO SCIACCA, AND JORGE P. ZUBELLI 

Abstract. We describe how the Harry Dym equation fits into the the bi-Hamiltonian 
formalism for the Korteweg-de Vries equation and other soliton equations. This is 
achieved by means of a certain Poisson pencil constructed from two compatible 
Poisson structures. We obtain an analogue of the Kadomtsev-Petviashivili hierar¬ 
chy whose reduction leads to the Harry Dym hierarchy. We call such system the 
HD-KP hierarchy. Then, we construct an infinite system of ordinary differential 
equations (in infinitely many variables) that is equivalent to the HD-KP hierarchy. 

Its role is analogous to the one played by the Central System for the Kadomtsev- 
Petviashivili hierarchy. 


1. Introduction 

The ubiquitous presence of completely integrable models, such as the nonlinear- 
Schrodinger equation, the Korteweg-de Vries equation, and the A-wave system, in a 
wide variety of physical phenomena explains part of the attention soliton equations 
have received during the last thirty years. Another reason for such attention is 
the multitude of deep connections of soliton equations with different branches of 
mathematics. They range from algebraic geometry to numerical analysis and include 
operator theory as well as symplectic geometry. 

Directly connected with the Korteweg-de Vries (KdV) equation one can find an 
intriguing equation known as the Harry Dym (HD) equation 

(1-1) qt = 2(l/\/(1 -|- q))xxx 

or equivalently 

(1 • 2) Pt P Pxxx 

after the substitution p = —(1 -|- 

Equation (p..l|) was discovered in an unpublished work by Harry Dym |I[|, and 
appeared in a more general form in works of P. C. Sabatier [^, ^, ^. It possesses 
many of the typical properties of soliton equations. See 0. B 0. B 0 and refer¬ 
ences therein. More recently, its relations with the Kadomtsev-Petviashvili (KP) and 
modi£ed-KP hierarchy have been studied in detail by Oevel and Carillo 0- 
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In the present work we discuss the HD equation from the bi-Hamiltonian point 
of view and show that it is amenable to the systematic treatment developed in 
O, More precisely, we start from a pair of compatible Poisson structures, 


12 


which also play an important role in the KdV hierarchy, and construct a Casimir 
for the corresponding Poisson pencil. By virtue of the standard construction the full 
hierarchy of Hamiltonian flows commuting with HD emerges. We then construct, in 
an elementary way and without resorting to pseudo-differential operators, a hierarchy 
of commuting flows which by reduction leads to the HD hierarchy. This reduction 
process is totally analogous to the one that leads from the KP hierarchy to the KdV 
hierarchy, which justifies our calling such hierarchy the HD-KP hierarchy []. Finally, 
we write down a system of infinitely many ordinary differential equations that admits 
the HD-KP hierarchy as a suitable reduction. 

The plan for the present work is the following: 

In Section o we review the general definitions of Poisson geometry, the bi- 
Hamiltonian theory, and the construction of the KdV and KP hierarchies from this 
view point. This will serve as a blue-print for the HD theory developed later in this 
work and will also help us emphasize the parallels between the HD-KP hierarchy and 
the usual KP hierarchy. 

Section is devoted to exhibiting the bi-Hamiltonian structure of the HD equation 
and to constructing an inhnite number of commuting flows from a solution of a Riccati 
equation. They will constitute the HD hierarchy. 

Section ^ describes the construction of the HD-KP hierarchy, i.e., a hierarchy that 
plays the same role for HD as the KP plays for the KdV. In our bi-Hamiltonian 
formalism, this hierarchy is depicted as an infinite number of conservation laws. The 
currents for such conservation laws are studied in Section In that section we also 
show that we can use such currents to fully describe the evolution of the quantities 
in the HD-KP hierarchy as ordinary differential equations. This leads to the notion 
of the Central System, which was also successfully used in earlier works dealing with 
the KdV and hidden-KdV hierarchies |^, |T^, pT| . 


1.1. Preliminaries: The bi-Hamiltonian Approach to Soliton Equations. In 

order to make the present work as much self-contained as possible we shall present 
in this sub-section an overview of the bi-Hamiltonian approach to the KdV and 
KP hierarchies. For a more detailed view and motivations the reader is referred to 


n. Hi lii Ii5], 0,0, [H, lii 


Our approach is based on the concept of Poisson manifolds, which is the natural 
setting to study Hamiltonian systems. 


Definition 1.1. Let M. he a differentiable manifold. A Poisson bracket on M. is an 
anti-symmetric "^-bilinear map 

(1.3) X ^ 

with the following properties: 

^In fact, Carillo and Oevel construct in [0 such HD-KP hierarchy making use of pseudo¬ 
differential operators. 
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1. Jacobi identity: {{F, G},H} + {{H, F},G} + {{G, H}, F} = 0. 

2. Leibniz rule 

{FG, H} = F{G, H} + {F, H}G. 

A Poisson manifold is a differentiable manifold endowed with a Poisson bracket. 

Starting from a Poisson bracket one can introduce a tensor P of type (2,0), which 
we shall consider as a map from T*Af to TfA, dehned by 

(1.4) {dG,PdF) = {F,G}. 

This tensor is called the Poisson tensor associated with {■,■}. In a Poisson manifold, 
the vector held Xh = {H, ■} = PdFd is called the Hamiltonian vector field associated 
with H. The map H i—> Xh, assigning to a function H its Hamiltonian vector held 
Xh, is a Lie algebra homomorphism: 

(1.5) X{f^g} = [Xf,Xg\. 

It follows that {F,G} = 0 implies [Xh,Xg] = 0; in other words, Hamiltonian vector 
helds associated with functions in involution commute. A Casimir function is a 
function H such that = 0, i.e., a function that is in involution with any other 
function on AI. In terms of the Poisson tensor, a Casimir function is a function 
whose diherential belongs to the kernel of P: PdH = 0. 

Definition 1.2. Let AA be a manifold endowed with two Poisson brackets, {•, - jo and 
- ji. These brackets are said to be compatible if any linear combination 

(1.6) {F,G'}a =^F,G}i-A{F,G}o 

is still a Poisson bracket. A bi-Hamiltonian manifold is a manifold M endowed with 
two compatible Poisson brackets. In this case the bracket {•, -Ia is called the Poisson 
pencil defined by {•, - jo and {•, .}i on AA. 

Since antisymmetry and Leibniz rule are automatically satished for {•, -Ia; the only 
remaining condition is the Jacobi identity. There are other equivalent forms for the 

def 

compatibility condition. For example, let us set for F G G°°{AA), Xp = P^dF and 

Yp = PidF. It is not hard to show that the compatibility condition can be rewritten 
as 

Lx,Pi + Ly,P, = 0 VF e G^{M), 
where Lx denotes the Lie derivative along the vector held X. 

Definition 1.3. Let {AA, Pq, Pi) be a bi-Hamiltonian manifold. A vector field X on 
AA is said to be a bi-Hamiltonian vector held if it is Hamiltonian with respect to both 
Pq and Pi, that is, if there exist two functions Hq and Hi such that 

(1.7) X = PodHi = PidHo. 

The next dehnition is the essential component behind Lenard’s construction. The 
latter is widely used to show that certain Hamiltonian models, such as the KdV, 
possess inhnitely many conserved quantities that are in involution with one another. 
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Definition 1.4. Let {Ai, Pq, Pi) be a hi-Hamiltonian manifold. A hi-Hamiltonian 
hierarchy is a sequence {Hk}k>o of functions on A4 fulfilling Lenard’s recursion re¬ 
lations 

( 1 - 8 ) {■, Hk+i}o = {■, Hk}i, k>0, 

and the additional condition = 0 . In terms of Poisson tensors: P^dH^+i = 

PidHk, with Hq a Casimir function of Pq: P^dHn = 0 . 

A bi-Hamiltonian hierarchy immediately gives rise to an inhnite sequence of bi- 
Hamiltonian vector helds 


( 1 . 9 ) Xk — PodHk — PidHk-i. 

They all commute, since the functions of a bi-Hamiltonian hierarchy are in involution 
with respect to both Poisson brackets. 

Remark 1.5. Let {Hk} be a bi-Hamiltonian hierarchy. Let us introduce the function 

( 1 . 10 ) H{\) = Y,Hk\-’^; 

k>0 

then, as an immediate consequence of the dehnition of bi-Hamiltonian hierarchy, we 
have that H{X) is a Casimir function of the Poisson pencil {•,-Ia- Vice-versa, let 
H{X) be a Casimir function of {•, -Ia which can be developed in a Laurent series: 
H{X) = J2k>-n^kX~^, with a suitable n; then the coefficients {Hk} form a bi- 
Hamiltonian hierarchy. 

From the above simple remark we get a method for seeking bi-Hamiltonian hierar¬ 
chies, namely: seek the Casimir functions of the Poisson pencil, which are deforma¬ 
tions of Casimir functions of Pq. This method has displayed its effectiveness in the 
case of a number of integrable models, such as the KdV hierarchy, the Toda lattice 
hierarchy, etc. See [^, 0, [T^, |^ . 

Remark 1.6. Let {Hk} be a bi-Hamiltonian hierarchy, and Xk = PodHk = PidHk-i 
the corresponding bi-Hamiltonian vector fields; let us introduce the functions 

k 

(1.11) R'=(A) 

j=0 

Then it is easy to show that, for all A, the vector fields Xk are Hamiltonian with 
respect to Px, with Hamiltonian H^~^{X): 

( 1 . 12 ) Xk = PxdH’^-\X). 

Therefore, the vector fields Xk are not only Hamiltonian, but they are Hamiltonian 
with respect to any bracket of the Poisson pencil. 


Let us briefly recall now the construction of the KdV and KP hierarchies of the 
bi-Hamiltonian point of view. For more details the reader is referred to |0, IT^ ITB 
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Let U be the infinite dimensional manifold of scalar valued C'^-functions on 
endowed with the Poisson brackets 

(1.13) {F, G}o = 2 / dF{dG)xdx 

Js^ 

(1.14) {F,G}i = ^ [ {dF{dG)xxx-4:u{dG)x-2uxdG)dx, 

2 J51 

so that the Poisson pencil is 

(1.15) (^Fx'ju'v ~^Vxxx F 2(u F ■^)vx F Ux'^- 

In these formulas, u G W is the point, and n is a co-vector at m, the pairing between 
V and a generic tangent vector ii being given by 


(1.16) {v,u) = / v{x)u{x)dx. 

Js^ 

In order to show that the Poisson pencil (|1.15|) admits a Casimir function we have to 
find an exact 1-form n(A) = dH{X) on U, such that 


(1.17) 


Pxv{X) = ~v{X)xxx F 2{u X)v{X)x F Uxv{X) = 0. 


Let v{X) be the unique solution of the equation 
(1.18) 


- ^'VVxx + (m + X)v‘^ = X 


of the form n(A) = 1-1- UjA L Then it can be shown that n(A) is an exact 
1-form. Indeed, if we set X = z‘^ and 


(1.19) 

so that 

( 1 . 20 ) 


V 2v 


hx F — u F y 


we have that the 1-form v is the differential of the function 


( 1 . 21 ) 


H = 2z 


h dx 


'Si 


which, therefore, is a Casimir of the Poisson pencil. Thus the Casimir functions of 
the Poisson pencil ( |1.15| ) can be computed by solving the Riccati equation ( p..20D . 

The basic object in the previous discussion is the map J from the manifold U to the 
space of Laurent series assigning to the point u the corresponding Hamiltonian 
density h{u), i.e., the solution of the Riccati equation ( |1.2CI|) . Since J allows us 
to compute the Hamiltonians, we call it the momentum map associated with the 
Abelian algebra acting on U by means of the KdV hierarchy. As a consequence of 
the involutivity relations between the functions Hj, when u evolves according to the 
KdV flows h{u) must obey conservation laws of the form 

dh _ 
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where the are suitable current densities. The important point is that one can 
write the directly in terms of a generic Laurent series of the form 

(1.23) h{z) = z + hjZ~K 

J>i 

Then the equations ( p..22|) become an infinite systems of partial differential equa¬ 
tions in the which is equivalent to the KP hierarchies. These ideas will be 

implemented in detail for the case of the HD equation in the next section. 


2. The Harry Dym Hierarchy and its bi-Hamiltonian Structure 

In this paragraph we start studying the HD hierarchy from the bi-Hamiltonian 
view point in analogy with the scheme developed earlier for KdV equation. Again, 
let lA be the space of scalar valued C'°°-functions on endowed with the Poisson 
brackets lEl, ^ 


( 2 , 1 ) 

( 2 . 2 ) 


dF{dG)xxx dx 


{F, Gji = — / {2u{dG)x + UxdG) dx. 
Js^ 


The corresponding Poisson tensors are given by 


(2.3) 

(2.4) 

SO that the Poisson pencil is 


‘^^dx Ux^ 


( 2 . 6 ) 


(^F\^uV ^xxx T •^(2'UUa; -j- Ux'v'). 


It is straight-forward to show that P\ is a Poisson tensor, so that the above defined 
brackets are compatible and W is a bi-Hamiltonian manifold. 

Our immediate goal is to show that, as in the KdV case, the Poisson pencil (p.5|) 
admits a Casimir function in the form of a Laurent series. This allows us to write 
a hierarchy of integrable vector fields for which HD is one of them. The important 


point is that the search for a Casimir function for (|2.5|) is also related to solving a 
Riccati equation. To wit, we start by finding an exact 1-form w{\) = dK{\) on U, 
such that 


( 2 . 6 ) 


Paw(A) ^ + X{2uyj{\)^ + u^'u-iX) ^ 0. 


We integrate once this equation noticing that 
(2.7) W (^^Wxxx + \{2uwx + ^ 

and therefore we consider the equation 


-^wl - -wwxx + 


( 2 . 8 ) 


-7^1 - i^wWxx + = A . 

4 2 
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This equation admits a unique solution of the form w{X) = Wo + X^j>i WiX \ and the 
coefficients Wi can be computed iteratively. In order to show that w{X) is an exact 
1-form, we set X = and remark that equation ( |2.1j| ) can be written in the form of 
a Riccati equation 


(2.9) 

For this reason, we set 

( 2 . 10 ) 

so that (|2.9| ) takes the form 

( 2 . 11 ) 


2) Wx\ ( Z Wx\^ 9 

W 2w/x \w Zw/ 


k[z) = - + —, 
w 2w 


kx + K= uz^. 


Then we consider a curve u(t) in W, and denote by w(t) the solution of (pf8|) at the 
points of this curve. By differentiating equation (|2.11|) with respect to a time t, we 
get z'^u = kx + Zkk and therefore 


( 2 . 12 ) 

{w,u) = 


Js^ 


w{x){kx + Zkk) dx = — [ { — -Wx + kw)kdx = 


Js^ 


dt Jg! 


k dx 


where we used equation (|2.10|) in the form 
(2.13) 


-Wx + kw = z. 
2 


Equation (|2.12|) shows that the 1-form w is the differential of the function 

(2.14) K = - [ kdx 

z Jsi 

which, therefore, is a Casimir of the Poisson pencil. 

Thus the Casimir functions of the Poisson pencil ( |2.5| ) can be computed by solving 
the Riccati equation ( |2.11|) ; and note that its unique solution has the asymptotic 
expansion in powers of z. 


(2.15) 


k{z) = k^iz + kj 


y-3 


j>0 


It can be shown that k 2 j is a total derivative for all j, so that K is a Laurent series 
in A. For their importance, we summarize these facts in the following 


Proposition 2.1. Let k he the unique solution of the Riccati equation \2.11\ ) 


kx k — uz 


admitting the asymptotic expansion (\2.1^ . Then, its integral K{X) given by (\2.14\ ) 
is a Casimir function of the Poisson pencil (ig. 
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Now we have an algebraic technique to compute the Casimir function. Before 
proceeding, we give an example how it works for the hrst orders of 2 ;. First, write 
the Riccati equation (2.11) 


( 2 . 16 ) 


k'^i = u 

(z) k-i^x + 2kok-i = 0 

( 1 ) ko^x kQ 2k-iki = 0 

(^z ki x R 2/cq/ci -|- 2k—\k2 — 0 

k2,x + kl + 2kok2 + 2k-ikz = 0 


for the hrst order coefficients. Thus we can solve (|2.16|) recursively, 

(2.17) k_i 

( 2 . 18 ) fco 

(2.19) ki 

( 2 . 20 ) 

At this point we have the expression in u and its derivatives of the coefficients kj 
of k{z). Following the previous proposition, we can construct the Casimir K{\) = 
2z~^ k dx of Pa- Hence, K = Kq + 7^2+ ..., whenever 

Kj = 2 [ kj-i . 

Js^ 

The functions K 2 j form a bi-Hamiltonian hierarchy and the vector helds of the hier¬ 
archy are now dehned by X 2 j+i = PodK 2 j = PidK 2 j- 2 - A simple calculation shows 
that the vector held corresponding to the Harry Dym equation (|1 . 1|) after the change 
u = 1 + q and a simple time rescaling is precisely 

^3(“) = ■ 

Finally, we conclude this section observing that all the HD equations are Hamil¬ 
tonian not only with respect to both Poisson tensors Pq and Pi, but also with respect 
to the Poisson pencil P\. Indeed, from a simple computation (see also Remark PD, 
we have that 

( 2 . 21 ) X 2 J +1 = PxdK^^^-^\X) , 

where 07^(21-2) 


u 

1 Ux 

4 u 


uz 


32 K 


1 

+ X 




8 uPu 


^Here we denote by the (S'(z))+ the non-negative powers of a Laurent expansion S{z) in z. 
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3. The Harry Dym-KP hierarchy 


The flows deflned by equation (|2.21|) are naturally plunged into a hierarchy of 
commuting flows with an infinite number of fields in a similar way that the KdV is 
plunged into the KP hierarchy. In fact, Carillo and Oevel have obtained such 
Harry Dym-KP hierarchy with the help of pseudo-differential operators. We close 
this section explaining how to relate the Harry Dym-KP hierarchy constructed in 
rT| with the one constructed here following the same techniques explained in . 


One of the key remarks of the present note is that we can construct such hierarchy 
“from first principles” as the one that governs the evolution of k{u), solution of (^.111) , 
as u evolves according to the Harry Dym hierarchy. The basic object is once again 
the map from the manifold U to the algebra of Laurent series, deflned by the Riccati 
equation ( p.ll|) . Since a solution k{z) of it allows us to compute the Hamiltonians, 
we call it the momentum map associated with the Abelian algebra acting on U by 
means of the HD hierarchy. As a consequence of the involutivity relations between the 
functions Kj, when u evolves according to the HD flows, k{u) must obey conservation 
laws of the form 

dk 


dt 




2j+l 


for some suitable current densities Our problem is to And a nice expression 

for the currents which can be computed following the ideas developed in 
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for KdV. From equation (|2.11|) we have that 


(3.1) 


{dx + 2k) 


dk 

dti 


= z 


du 

dtj' 


while the HD equations ( |2.21|) can be written as 
(3.2) 


dt 


3ij 1 

. = [dx + 2k){-dx){-dx + 2k){X=w) + , 


2j+l 


Since 

(3.3) 

Therefore 


^Wxx. 


+ XuWx + UxW-{dx + 2k){^dx){ dx + 2k)- 


V. 


{dx + 2k) = z\dx + 2k){\dx) ^w)+^ + 2k{z^^ ^w)+) 

= (dx + 2k)dxZ^ + hiz^^-^v)+) , 


which implies 

(3.5) = -^{X^v)+,x + h(A%)+. 

This is the formula for the current densities that appear in our local conservation 
laws. Observe that we can put since the even times are trivial, and we 

extend the definition of the currents for every j. 
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The fundamental observation allowing to generalize this procedure to the HD-KP 
theory is that we can compute the currents directly from k, without computing 


the series w{X), in perfect analogy with the KdV theory 


A key object in our construction is the notion of Fad di Bruno polynomials. They 
are defined as + ky ■ 1, for j > 0. 

Proposition 3.1. The above defined currents fori > 2, are determined directly 
from k, and in particular they satisfy 

1 . = z^ + 0{z) 

2. G K'f^, where K'f^ =‘^ ■ ■ ■) (linear span over -functions). 


Sketch of the proof: For I even, there is nothing to show. To show the hrst item 
above for / = 2j + 1 we write {z‘^^~‘^w)_^_ = z‘^^~‘^w — {z^^~‘^w)_ and make use of 
the relation —\wx + kw = z. As for the second item, we recall that {z'^^~^w)^ = 
and thus 

k=0 k=0 

To conclude the proof of the second item, it remains to show that z'^^k^^^ G K'(^ and 
hence 

z^k^i) = z^dx + kYil) = {dx + k)\z^) G Kl 
since {dx + k)K, c K'f. Thus, for example. 


U U 


+ ’ 


and so on. Q. E. D. 

The previous proposition is the central point to formulate the KP theory for the 
HD equations. It means that we can compute the currents with j > 2, starting 
with any Laurent series k{z) of the form 

(3.6) k = k-iz + kiz~’‘ 

l>0 

by following the two requirements in items (1) and (2) of Proposition H and not 
necessarily those satisfying ( p.ll| ). For example, it is possible to write 

(3.7) = cik^^^ + • • • + , 

with Cj coefficients, that depend on the coefficients of A;, determined by the asymptotic 
expansion. At this point the local conservation laws became evolution equations 
in the space of Laurent series of the form (|3.6|) and live independently of the HD 
equations. 


Definition 3.2. The Harry Dym-Kadomtsev-Petviashvili (HD-KP) hierarchy is the 
system of vector fields acting on the space of Laurent series of the form defined 
by 


dk 

dti 


dxK^''> , 1 = 2,3,-■■ 
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where the conserved densities are defined by the prescription of items (1) and 
(2) of Proposition 


To obtain the HD hierarchy dehned above all that we need to do is to impose the 
constraint 


(3.8) 




Indeed, it is easy to see that (p.8|) is equivalent to the Riccati equation for Harry Dym: 
express, according to (|3.7|) , the current in the basis of Faa di Bruno polynomials, 
i.e., 

R:(2) = ; 

then (^) becomes the Riccati equation ( 2.11|) if we put u = k-i^. This implies that 
_ ^ 2 k k = 1,2, ■ ■ ■, and thus all the even flows are stationary. 

In O a hierarchy of commuting flows was constructed, making use of the no¬ 
tion of pseudo-differential operators, that includes the Harry Dym flow as one of its 
reductions. To do that, they consider the operator 


(3,9) 


def 


^Dym 


— Q-ldx + Q'O + + ■ ■ ■ 


The hierarchy is then defined in Lax form as 

^.^Dym 


(3.10) 


dti 


(-^Dym) >2 ’ -^Dym 


,J>2 


where (S ')>2 means the projection on the differential operators of order > 2 of the 
pseudo-differential operator S. 

This definition is equivalent to ours, after a change of variable which follows the 
T3[|. Start by defining the negative Faa di Bruno polynomials with j < 0, 


one m 


which can be obtained from k^^'^ = 1 by solving backwards the equations 

kU+i) = kij) + 

and expand 2 ; on the basis 

2 ; = Q-ik + qjk^~^\ 

j>0 

This establishes a change of variables between {kj}j>-i and {qj}j>-i- Then it can be 
shown that equations ( p.lOl) coincide with the HD-KP hierarchy given in Definition 


3.2 


4. The Evolution of the Currents 

We shall now study in more detail the evolution of the currents when k evolves 
according to the HD-KP hierarchy. Our hrst remark is that we can interpret the HD- 
KP equations as the commutativity conditions 

[dt^+K^^\dx + k] =0 
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of the operators dt^ + and dx + k. Since G K'^ and K'^ is invariant with 
respect to dx + k, it follows that 

(4.1) {dt^ + = {dt^ + {dx + kyi = {dx + ky e <, 

for every I > 0 and j > 2. Therefore, 

(4.2) {dt^ + {K^^^) e K'l . 

As in the KdV case, we observe that the currents Kd"> form a basis of K'^. This fact 
allows us to look at the equations (|4.2| ) as a combination of the currents with suitable 
coefficients. In more detail, if we write 

(4.3) = z* + K\z + Kiz-^ , i > 2 , 

l>0 


then, equations (f4.2|) can be written in the explicit form 
dK^i) 


dtj 




i-2 


J-2 


(4,4) +5^A7ii''-'> + 5; forz,J>2, 

1=1 i=i 

which gives the evolution of the currents. A similar expression can be found for the 
evolution of k. Indeed, ( [4.1|) entails that 


so that 


dk 

dtj 


+ kK^^^ = + koK^^'> 


i-2 


(4.5) 


+ Y kiK^^-^\ for j > 2. 


1=1 


Of course, this is not different from the dehnition of the HD-KP hierarchy, but it 
does not explicitly contain x-derivatives. Furthermore, the following equations also 
hold 


dxK^^^ + kK^^^ 


(4.6) 

(4.7) dxk + k'^ 


i-2 

+Ki^k.iK^'^'^ + Y for j > 2, 

1=1 


coming from the fact that the operator dx + k sends the span K'j^ 

K- 


{kd^)j>i into 
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In the next section we shall see that, in analogy with the Central-System equations 
for the KP hierarchy, equations ( |4.4|) , ( |4.6| ), ([4.5|), and (|4.7|) will give rise to the 
Central-System equations for the Harry Dym hierarchy. 

We stress that the above equations underscore the fact that the evolu¬ 

tion of the Laurent coefficients of the currents and k undergoes a dynamical 
system which does not involve spatial derivatives. We also remark that such system 
is different from the one obtained in since here k is not monic anymore. 


5. The Central System of the Harry Dym Hierarchy 


In the previous section we analyzed the evolution of the currents when the HD- 
KP equations hold. We notice that all the coefficients in the equations (|4.4H4.7|) are 
dependent on the coefficients kj of k through the Faa di Bruno polynomials as in 
Proposition |3.1| . The main point to pass to the Central System is to consider the 
currents (|4.3|) as arbitrary elements of K'f independent of k. This is one step further 
to considering k independent of the Riccati equation when we passed from the local 
conservation laws to the HD-KP equations. In this new picture, the series k does not 
play any special role, and we rename it by 

+ 

l>0 


We have shown that the evolution of the currents is given by the equations 

(5.1) {dt^ + iP(') e kI , for l,j > 1. 

It is not hard to see that they can be rewritten in the form 

(5.2) , for> 1, 

where we denote by and the projections of the space L of all Laurent series 
onto the subspaces K”^ and K'f respectively, defined by the following decomposition 

(5.3) L = z^-^\ ...) = kI®k1. 


We call ( p^ the HD-Central System (HD-CS) associated to the HD-KP theory. We 
conclude the present section by analyzing a few of its properties. 

The first remark is that the evolution of the currents given by the Central-System 
equations for the Harry Dym hierarchy form a system of nonlinear ordinary differ¬ 
ential equations of Riccati type in the components of the fields for i > 1. The 
second remark is that the flows defined above commute with one another. To show 
that, we start with the remark that along the HD-CS flows (|5.2| ), we have that 


dti dtj 

Then, we compute the commutator of two vector fields of the HD-CS hierarchy on a 
generic current Let us call such fields say Xj and X;, and notice that 


[X„X,] (x«) 


d d 

dtj dti dti dtj 


'A 

9tj 


+ , 


A 

dti 


+ 7^(0 
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Because of the specific form of the Laurent series of the above commutator 
belongs to On the other hand, because of the invariance condition ( b.2|) that 
dehnes the flows, we have that the commutator belongs to iL", and hence vanishes. 

As a consequence we can apply the spatialization process, already described in 
earlier works 


2^. This procedure is based on the following: 


Remark 5.1. Assume we have a family {X^} o/commuting vector fields on a mani¬ 
fold. Then, two reductions processes are natural. The first one is the projection along 
the integral curves of one of the vector fields. The second one, is the restriction to 
the set of zeroes of a fixed vector field. 


It easy to see that, in particular, by setting ti = x and using the first flow to write 

all the equations in terms of k and its derivatives with respect to x we obtain 

the HD-KP hierarchy. This corresponds to the hrst reduction mentioned above. The 
second step is to consider the stationary flows with respect to the t 2 flow given by 
the constraint = z^. This gives back the Riccati equation and the HD hierarchy 
we saw in Section |^. 


6. Concluding Remarks 

In the present work we have carried out for the Harry Dym equation, the first 


steps of the program that was developed in [13, 18] for the KdV and KP hierarchies. 


Namely, we constructed the HD hierarchy, the HD-KP hierarchy and the correspond¬ 
ing Central System. The program in ||I3| , p!8| lead to the linearization of the flows 
in the KP hierarchy and the construction of polynomial r-functions This 

was achieved by using the notion of Darboux coverings and introducing the modified 
Central System. This would be a natural follow up of the present work. 
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